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ABSTRACT 


■vl 

GOES-VAS  satellite  retrievals  created  by  a  new 
Hsimul taneous1'  algorithm  are  examined,  for  a  case  in  July 
1982.  *>The  primary  goals  are  to  determine  the  spatial  and 
temporal  character i sti cs  of  the  data,  as  well  as  the 
magnitudes  of  mean  and  RMS  errors.  several 
integrated  and  nan-integrated  stability  parameters  are 
examined  to  assess  VAS's  capability  for  providing  useful 
information  in  convective  storms  forecasting. 

Results  indicate  that  the  new  algorithm  provides  a 
much  better  estimate  of  moisture  than  the  older  "physical" 
variety.  On  the  other  hand,  temperature  showed  no 
i mprovement . 

„  The  improved  moisture  information  has  a  positive 
impact  on  various  stability  indices.  A  VAS-derived  Lifted 
Index  appears  to  be  a  promising  convection  forecast  tool. 
Many  small  scale  features  in  satellite  imagery  were 
related  to  the  Lifted  Index  analyses.  Also,  the 
VAS-derived  values  provided  information  that  could  not  be 
obtained  from  conventional  surface  data.  The  K  Index  and 
Total  Totals  Index  also  provided  forecast  data;  however, 
their  resolution  did  not  appear  as  good  as  those  of  the 
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1.  Introduction 


With  the  launch  of  the  -first  Visible  In-frared  Spin-Scan 
Radiometer  (VISSR)  Atmospheric  Sounder  (VAS)  onboard  a 
geostationary  platform  in  the  fall  of  1980,  forecasters 
have  had  a  significant  increase  in  the  amount  of  data 
available  for  the  analysis  and  prognostication  of  both 
large  and  small  scale  weather  phenomena  such  as  severe 
storms  and  mid-latitude  cyclones.  The  12  VAS  channels 
measure  infrared  radiation  in  bands  from  3.5  to  15  pm. 
Their  information  can  be  used  in  a  variety  of  ways.  For 
example,  channel  imagery  can  be  used  independently,  such 
as  done  with  the  6.7  pm  water  vapor  band  which  is  used  to 
observe  mi d-tropospher i c  humidity  (e.g.,  Petersen  ejt  al .  . 
1982).  Other  applications  of  6.7  pm  imagery  include  the 
generation  of  horizontally  enhanced  fields  of  middle  level 
mixing  ratios  (e.g.,  Steranka  aj_.  ,  1973),  wind 
determination  by  tracking  vapor  features  (e.g.,  Stewart  et_ 
al . .  1985),  and  inferring  the  character i sti cs  of  upper 
level  jet  streams  and  jet  streaks  (e.g.,  Rodgers  al .  . 
1976;  Raymond  et^  al_.  ,  1981;  Stewart  and  Fuelberg,  1986). 
For  deriving  vapor  content  nearer  the  surface,  Chesters  et_ 
al .  (1983)  used  the  "split  window"  technique,  a 
combination  of  VAS  channels  7  and  8.  When  employed  with 


channel  in-formation  describing  other  levels  of  the 
atmosphere,  areas  of  convective  instability  (dry  mid-level 
air  over  moist  low-level  air)  could  be  calculated 
(Petersen  et_  al .  ,  1984)  . 

If  measurements  from  all  VAS  channels  are  utilized,  the 
radiances  can  be  transformed  into  vertical  temperature  and 
moisture  profiles  based  on  knowledge  of  the  associated 
vertical  weighting  functions.  There  are  currently  two 
different  approaches  for  creating  these  soundings.  One 
method  requires  the  development  of  statistical 
relationships  between  measured  channel  radiances  at 
radiosonde  locations  and  the  actual  radiosonde 
measurements.  These  regression  equations  then  are  used  to 
calculate  atmospheric  parameters  at  the  various  levels  of 
different  sites.  Conventional  surface  data  also  are 
frequently  used  in  the  equations  (Lee  et_  al_.  ,  1983)  . 
Studies  using  the  "statistical"  retrievals  demonstrated 
reasonable  agreement  with  soundings  from  a  special 
mesoscale  radiosonde  network,  especially  in  the  lower 
troposphere  (Chesters  et  ai_-  ,  1984).  The  ability  of 
statistical  retrievals  to  describe  the  atmosphere's 
current  moisture  content  has  been  examined  extensively. 
For  example,  investigators  have  noted  that  they  can 
provide  valuable  information  about  the  pre-convecti ve 
environment  (Petersen  et^  ai_.  ,  1984).  Petersen  al  - 
(1982)  demonstrated  that  VAS  stability  calculations 


provided  valuable  in-formation  between  radiosonde  data 
points,  whereas  Mostek  et_  al  -  (1986)  noted  that 
regression-deri ved  stability  parameters  exhibited 
accurate  mesoscale  resolution  and  continuity. 

The  second  method  of  retrieval  calculation  involves 
inverting  the  radiative  transfer  equation.  This  "physical" 
algorithm  (Smith,  1983)  uses  an  iterative  convergence 
scheme  to  determine  the  sounding  profile  which  fits  the 
measured  channel  radiances  with  the  least  amount  of 
residual  error.  The  iterations  loop  between  the 
temperature  and  dewpoint  profile,  making  each  dependent 
upon  the  other.  A  study  of  mesoscale  air  masses  based  on 
these  physical  retrievals  (Zehr,  1985)  suggested  that  VAS 
information  could  provide  mesoscale  convection  forecasts. 
Leftwich  (1985)  found  hail  forecasts  using  VAS  profiles 
were  timelier  and  more  accurate  than  radiosonde  generated 
profiles.  Operational  use  of  physical  retrievals  at  the 
National  Severe  Storms  Forecast  Center  (NSSFC)  has  been 
documented  by  Anthony  and  Wade  (1983).  One  unfortunate 
aspect  of  VAS  physical  profiles  is  a  large  2  -  4  “c  diurnal 
temperature  variation  extending  well  above  the  boundary 
layer  (Fuel berg  and  Funk,  1987). 

Investigations  into  the  accuracy  of  the  two  VAS 
retrieval  methods  have  been  conducted  by  several 
researchers.  Comparisons  with  radiosonde-derived  "ground 
truth"  data  (Jedlovec,  1985)  have  indicated  bias  errors 


less  than  4  C  -for  temperature  and  3-6.5  C  -for  dewpoint 
profiles  in  the  case  of  physical  retrievals.  Errors  less 
than  3.5*C  (temperature)  and  less  than  2  °C  (dewpoint) 
were  associated  with  the  regression  algorithm.  Physical 
temperature  retrievals-  had  slightly  smaller  Root  Mean 
Square  (RMS)  errors  (1.7*0  than  the  regression  retrievals 
(2.0*0.  In  contrast,  however,  the  regression  dewpoint 
profiles  had  smaller  RMS  errors  (4.8*0  than  the  physical 
retrievals  (5.9*0.  Results  of  both  algorithms  also  have 
been  evaluated  using  statistical  structure  functions 
(Fuelberg  and  Meyer,  1985).  Random  errors  in  the  VAS 
profiles  were  comparable  to  those  from  radiosondes,  with 
the  physical  algorithm  showing  slightly  better  results. 
Similarly,  horizontal  gradients  from  the  physically 
derived  soundings  were  more  meaningful  than  those  from  the 
statistical  approach. 

A  new  physical  retrieval  algorithm  (Smith  et  al_.  , 
1986;  Smith  and  Woolf,  1984)  has  been  devised  to  correct 
some  of  the  documented  problems  of  the  original  physical 
formulation  (Smith,  1983),  especially  the  large  errors  in 
moisture  content.  This  new  method  uses  a  simultaneous 
solution  system  rather  than  an  iterative  scheme. 
Temperature  and  moisture  calculations  are  made 
simultaneously,  rather  than  being  dependent  upon  each 
other.  Thus  far,  only  a  few  studies  have  utilised 
soundings  from  this  recent  approach.  Wade  at  al .  (1985) 


presented  several  thermodynamic  parameters  generated  in 
near  real  time  -for  the  NSSFC  during  May  1985.  Kitzmiller 
(1986)  tested  well  known  stability  indexes  generated  -from 
VAS  retrievals  and  found  some  predictor  combinations  to  be 
promising  for  severe  storm  forecasting. 

The  purpose  of  this  paper  is  to  examine  temperature  and 
dewpoint  profiles  created  by  the  new  simultaneous 
algorithm  for  a  period  in  July  1982.  Although  this  case 
was  previously  documented  by  Fuelberg  and  Funk  (1987) 
using  the  Smith  (1983)  physical  retrieval  algorithm,  the 
current  emphasis  differs  greatly  from  theirs. 
Specifically,  the  accuracy  of  temperature  and  dewpoint 
data  will  be  statistically  compared  to  radiosonde 
observations  and  to  the  older  physical  values  used  by  Funk 
(1985).  In  addition,  spatial  and  time  continuities  of  the 
data  fields  generated  by  VAS  will  be  examined.  After 
discussing  the  apparent  strengths  and  weaknesses  of  the 
new  algorithm,  numerous  traditional  parameters  for 
thunderstorm  forecasting  will  be  examined.  That  section 
will  focus  on  the  following  questions.  Did  VAS  retrievals 
provide  useful  temporal  or  spatial  information  which  was 
not  available  from  the  radiosonde  observations  (RAOBs)? 
What  was  the  influence  of  the  hourly  surface  data  in 
generating  the  low  level  VAS  temperature  and  moisture 
profiles;  that  is,  did  VAS  provide  information  not 
contained  in  the  hourly  surface  data? 


It  is  important  for  the  user  of  VAS  data  to  understand 
strengths  and  weaknesses  of  the  current  satellite 
retrieval  system.  The  intent  of  this  paper  is  to  help  shed 
additional  light  on  those  products. 


2.  VAS  Instrument 


The  VAS  instrument  is  a  passive  radiometer  capable  o-f 
making  simultaneous  visible  and  infrared  images  of 
radiation  upwelling  from  clouds,  the  earth's  surface,  and 
it's  atmosphere.  VAS  is  carried  onboard  the  latest  series 
of  geostationary  spacecraft  (GOES).  Compared  to  polar 
orbiting  platforms,  the  geostationary  configuration 
permits  a  much  faster  sampling  rate  for  observing 
particular  surface  locations.  The  instrument  spins  at  100 
rpm  from  west  to  east  and  steps  a  scan  mirror  from 
(usually)  north  to  south.  The  reader  is  referred  to  Menzel 
et  al .  (19B1)  for  a  detailed  description  of  the  VAS 
instrument  and  its  operational  modes.  Briefly  stated, 
however,  VAS  consists  of  8  visible  channel  detectors  and  6 
thermal  detectors,  the  latter  measuring  infrared  radiation 
in  12  channels  with  central  wavelengths  ranging  from  3.7 
to  IS  jj m  (See  Table  1,  Fig.  1).  Fields  of  View  (FOV)  of 
the  visible  detectors  possess  a  1  km  resolution  at 
satellite  nadir,  while  the  thermal  detectors  have  a  7  or 
14  km  resolution  depending  on  the  type  being  employed.  A 
filter  wheel  placed  in  front  of  the  detectors  selects  the 
spectral  channel  being  used  at  the  given  moment. 

Seven  of  the  12  spectral  channels  (Table  1,  Fig.l)  are 
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Table  1.  VAS  instrument  characteri sti cs  (a-fter  Smith,  1983). 
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thermal  C0Z  channels  <1-6  and 
measurements  can  provide  vertical  temperature  pro-files.  On 
the  other  hand,  there  are  only  3  channels  <7,  9  and  10) 

which  yield  moisture  data.  Finally,  channels  8  and  12  give 
surface  or  cloud  top  information. 

Knowledge  of  the  absorption  characteristics  of  several 
atmospheric  constituents  <C0Z  ,  Nz ,  Nz0,  NO,  CH,^ ,  02 ,  and 

O)  gives  the  weighting  functions,  d77d  In  p  of  the  12 
spectral  channels  displayed  in  Fig.  1.  Although  each 
channel  detects  radiation  over  a  considerable  atmospheric 
depth,  peaks  for  the  moisture  channels  are  relatively 
sharp,  due  to  the  rapid  decrease  of  water  vapor  with 
altitude  in  the  lower  troposphere.  Even  though  each 
channel  peaks  at  a  different  level,  the  weighting 
functions  have  considerable  overlap,  contributing  to  VAS's 
poor  vertical  resolution. 

VAS  has  3  operating  modes,  the  Dwell  Sounding  (DS) 
mode,  the  Mul ti -Spectral  Imaging  (MSI)  mode,  and  the  Dwell 
Imaging  (DI)  mode.  Each  can  use  any  of  the  visible  and 
infrared  channels.  Capable  of  receiving  instructions  from 
the  ground,  the  instrument  can  scan  any  subset  of  the 
earth's  disk,  with  scanning  time  dependent  upon  mode 
selection  and  size  of  the  observed  area. 

The  DS  mode  was  designed  to  improve  signal  to  noise 
ratios  of  the  radiance  measurements  so  that  sounding 
profiles  can  be  retrieved.  This  is  achieved  by  keeping  the 
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scanning  mirror  fixed  on  a  given  line  -for  a  predetermined 
number  of  instrument  revolutions  far  each  spectral  channel 
programmed.  Thus,  the  multiple  passes  over  the  same 
locations  provide  an  averaging  process  that  partially 
filters  random  error  in  the  data.  The  number  o-f  spins  per 
line  varies  with  channel  and  is  called  the  spin  budget 
(Table  1).  In  addition,  a  spatial  averaging  o-f  several 
horizontal  fields  of  view  also  is  used  to  improve  signal 
to  noise  ratios.  Of  course,  this  decreases  the  horizontal 
sounding  resolution,  often  down  to  near  75  km.  Due  to 
multiple  passes  along  each  scanning  line,  operational  time 
constraints  do  not  permit  the  entire  disk  to  be  observed 
in  this  mode.  About  30  minutes  is  required  to  observe  20 
degrees  of  lati  Jde  (Petersen  et  al . .  1982). 

The  MSI  mode  permits  larger  areal  coverage  and  more 
rapid  rescanning  since  only  a  limited  number  of  channels 
is  sampled  during  a  particular  scan  pass,  i.e. ,  up  to  4 
spectral  bands,  always  including  the  visible  channel.  In 
the  MSI  mode  the  earth's  disk  can  be  observed  every  30 
minutes.  Sounding  retrievals  are  not  prepared  from  MSI 
data. 

The  Dwell  Imaging  (DI)  mode  was  created  as  a  compromise 
between  the  DS  and  MSI  configurations.  Specifically, 
measurements  are  taken  over  all  12  channels;  however,  the 
spin  budget  is  smaller  than  for  the  DS  mode.  Nonetheless, 
sounding  retrieval  still  is  performed  although  results  are 


considered  not  as  reliable  as  those  -from  the  expanded 
budget.  The  D1  mode  permits  coverage  of  the  North  American 
continent  at  30  minute  intervals  (Smith  al_.  ,  1982). 

During  severe  storm  situations,  VAS  measurements  are 
pre-empted  by  Rapid  Interval  Scan  Operations  (RISOP)  that 
provide  more  -frequent  visible  and  infrared  window  channel 
(8)  imagery.  Thus,  retrieval  products  often  are  not 
prepared  during  severe  weather  events.  This  has  limited 
many  investigators  to  exploring  VAS  capabilities  during 
pre-convecti ve  and  non-severe  weather  cases. 


3.  Physical  retrieval  algorithms 


The  physical  retrieval  algorithm  involves  the 
conversion  of  measured  spectral  radiances  into  vertical 
profiles  of  temperature  and  moisture  by  using  the 
radiative  transfer  equation.  More  specifically,  least 
squares  analytical  solution  to  the  equation  transforms  a 
smooth  "first  guess"  profile  into  an  enhanced  final 
version.  Since  a  detailed  discussion  of  this  algorithm  is 
presented  by  Smith  (19S3)  only  a  summary  follows. 

The  radiative  transfer  equation  is  of  the  form 


ft,  (1) 

where  R  is  defined  as  the  mean  spectral  radiance  for 
wavenumber  V  and  local  zenith  angle  9.  B  is  the  "Planck" 
radiance  associated  with  temperature  T,  71  is  the 
transmittance  of  *he  atmosphere  above  level  p,  and  the 
subscript  s  denotes  the  opaque  surface  of  the  radiating 
atmosphere,  either  ground  or  cloud.  Both  B  and  Tf  are 
weighted  with  respect  to  the  spectral  response  of  the 


particular  channel.  The  model  treats 


B 


and  as 


monochromatic  functions,  thereby  simplifying  the 
computations. 

Input  data  are  assumed  to  be  measurements  of  R(V,®),  a 
transmittance  model  providing  f  functions,  and  estimates 
of  ps .  A  "first  guess"  temperature  profile  is  required  and 
can  be  constructed  in  a  variety  of  ways.  For  example, 
climatology  tables  can  be  used  to  specify  a  standard 
atmosphere  for  a  given  latitude  and  season.  Alternatively, 
a  numerical  analysis  or  forecast  (12  h)  from  the  NMC 
models  can  be  employed,  or  previous  VAS  profiles  can  serve 
as  the  "guess"  for  the  next  profiles.  In  addition, 
observed  surface  data  can  (and  usually)  are  blended  with 
these  starter  profiles.  Since  the  surface  (cloud  or 
ground)  is  assumed  to  be  opaque,  a  subroutine  is  included 
to  insure  this  requirement. 

A  "guess"  mixing  ratio  profile  is  created  from  the 
first  guess  temperature  sounding  by  assuming  SOX  relative 
humidity  (RH)  and  then  using  the  derivation  of  Smith  and 
Zhou  (1982) 


RH  =  RH0  exp 


R(N),e)  -  r£v>, 

jPs  Q(V,P;©)ap  ' 

o 


(2) 


where  R(  ^  , ©  )  is  the  observed  radiance  correspond! ng  to 


the  true  distribution  of  precipitable  water,  R(v,  RH^,,  0  ) 
is  a  radiance  calculated  for  an  assumed  distribution  of 
relative  humidity  RH0,  and  Q(v,p,6)  is  a  moisture 
weighting  function.  Mean  relative  humidities  are 
calculated  for  three  layers  bounded  by  800,  600,  and  300 
mb,  i.e.,  the  approximate  levels  of  weighting  function 
peaks  for  the  three  VAS  water  vapor  channels  (7,  9,  10). 
If  available,  values  of  surface  relative  humidity  also  are 
utilized,  otherwise  the  channel  7  derived  RH  value  is 
used.  A  "smooth"  water  vapor  profile  is  achieved  by 
linearly  interpolating  between  these  calculated  values. 
This  step  assumes  a  constant  RH  within  the  layer 
contributing  radiation  to  each  water  vapor  channel. 

□nee  the  guess  profiles  are  prepared  as  described 
above,  transmittance  functions  are  calculated  for  the  VAS 
thermal  (C02)  channels  (1,  2,  3,  4,  5,  6,  1 1 ) .  Thermal 
channel  radiances  then  are  estimated  by  using  (1)  and  (2) 
from  the  "guess"  temperature  and  smooth  moisture  profiles. 
Following  Smith  (1983),  differences  between  the  measured 
brightness  temperature  T0  <\>)  and  T0N(\)),  the  brightness 
temperature  which  would  be  measured  if  the  true 
temperature  profile  T (p)  equalled  the  estimate  T^tp)  can 
be  expressed  as 


Tb(V)  -  TB"(N))  =  CT(p)-T"(p)  3WM(V,p,0)  dp/p 
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where  N  denotes  the  iteration  step  and  WM ( ^  ,p,9>  is  the 
temperature  pro-file  weighting  function.  After  normalizing 
the  weighting  function  with  respect  to  the  expected  error 
of  the  brightness  temperatures ,  W“(V,p,0>,  an  iterative 
solution  of  (3),  proposed  by  Smith  (1970)  is  used: 

T--Mp>  =  T^p)  +  y  W-  (  V,p  ,0)  CT>  (  V)  -Ta  N  (V)  3  f  (4) 

J  W(V,p,6) 

where  j  is  the  spectral  channel  index. 

After  the  initial  "guess"  profile  has  been  modified  by 
(4) ,  a  new  set  of  C03  channel  radiances  is  calculated  in  a 
similar  manner.  Then,  a  convergence  test  is  applied  to 
determine  whether  the  new  temperature  profile  is 
significantly  different,  i.e.,  the  change  in  brightness 
temperatures  is  compared  to  the  "expected"  error.  When  the 
change  is  less  than  an  empirically  derived  noise  factor, 
the  new  profile  is  determined  to  be  the  final  version.  If 
the  profile  fails  the  test,  additional  iterations  are 
performed  until  convergence  is  achieved. 

The  iterative  scheme  unfortunately  produces  a  smooth 
temperature  profile.  To  partially  rectify  this  problem,  an 
enhancement  procedure  involving  least  squares  solutions  is 
applied.  Specifically,  it  is  assumed  that  the  difference 
between  the  true  profile  T(p)  and  the  retrieved  profile 
T~(p)  can  be  approximated  by  line  segments  expressed  by 

dT2, 


three  unknown  quantities  dT , , 


and  dT  5  as  shown  i  n 


Fig.  2.  Similarly,  the  water  vapor  profile  is  enhanced  by 
using  a  line  segment  approximation  to  the  differences 
between  the  true  and  retrieved  profiles. 

Based  on  the  above  mentioned  procedure,  an  enhanced 
water  vapor  profile  is  constructed  and  used  in  the 
transmittance  model  to  calculate  new,  updated  CO^ 
transmi ttance  functions.  Then,  a  new  enhanced  temperature 
profile  is  calculated  as  in  Fig.  2.  Next  an  enhanced  water 
vapor  profile  is  computed  again.  Finally,  dewpoints  are 
obtained  from  the  final  moisture  values. 


b.  Simultaneous  algorithm 


The  "new"  physical  retrieval  algorithm  employs  a 
different  procedure  for  solving  the  radiative  transfer 
equation  (Smith  et_  al_.  ,  1986;  Smith  and  Woolf,  1984). 
Instead  of  an  iterative  solution  for  the  profiles,  a 
simultaneous  retrieval  of  temperature  and  moisture  is 
performed.  As  in  the  original  method  described  above,  a 
"first  guess"  profile  is  required,  and  options  available 
for  its  construction  are  similar.  Only  highlights  of  the 
algorithm  are  presented  below. 

In  perturbation  form,  the  radiative  transfer  equation 
can  be  expressed  as 


t  i(Ss), 
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(5) 


Illustration  of  the  line  segment  approx¬ 
imations  used  to  model  the  fine-scale 
structure  deviations  of  (a)  atmospheric 
temperature  ana  <b>  water  vapor  from  their 
vertically  smoothed  estimate  (after  Smith, 
1983). 


where  dependency  on  angle,  pressure,  and  -frequency  are 
eliminated  -for  simplicity.  The  perturbations  (£)  are  with 
respect  to  an  a-priori  condition  (e.g.,  a  "-first  guess" 
pro-file).  Smith  et  al_.  (19B6)  assumed  that  the 

transmittance  perturbation  is  dependent  only  on  the 
uncertainty  in  the  column  o-f  precipitable  water  (U)  . 
Integration  o-f  (5)  and  application  o-f  the  assumption 
produce: 


ST 


"Jsu  dp  Jo  “  J  3p(dB/dT*)Ap  +  ^-r»)VSi(6) 


where  T**  is  the  brightness  temperature,  T  is  temperature, 
and  Ts  is  the  surface  skin  temperature. 

To  solve  for  the  perturbations  from  a  set  of  spectrally 
independent  radiance  observations,  arbitrary  pressure 
functions  fd(p)  are  used  for  JTj  ,  S  U(p)  ,  and  5  T (p)  in  (6). 
Smith  al_.  <19B6)  selected  the  profile  weighting 
functions  (d27d  In  p)  of  the  radiative  transfer  equation 
as  pressure  basis  functions.  Substituting  these  functions 
into  (6)  yields  a  spectral  radiance  observation,  $ Tj  ,  for 


the  set  of  K  channels: 
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Other  quantities  needed  to  calculate  the  ±  matrix  are 
estimated  -from  the  -first  guess  sounding.  The  exception  is 
dtVdU  which  is  obtained  from  a  mixing  ratio  profile  with 
values  at  each  level  being  half  those  of  the  first  guess 
profile. 

In  matrix  form  (7)  becomes 


±*  =3*, 


where  t~  is  a  row  vector  of  K  brightness  temperature 
observations,  0(.  i  s  a  row  vector  of  M+l  coefficients,  and  4 
is  a  matrix  having  dimensions  <KxM+l>.  A  least  squares 
solution  of  <9)  gives 


Oi  =  (STI  +  XT)"'  2xt* 


where  (  ) T  denotes  matrix  transposition,  <  )_l  indicates 

matrix  inversion,  and  Y  is  a  scaler  (usually  0.1) 


multiplied  by  the  identity  matrix  I  and  added  to  the  fTf 
matrix  in  order  to  stabilize  the  matrix  inverse.  After  the 
ds  are  determined,  $TS,  £q,  and  ST  are  known  (by 
knowledge  of  the  original  basis  functions).  These  values 
then  are  added  to  the  first  guess  sounding  to  achieve  a 
final  profile. 

There  are  two  passes  through  the  retrieval  algorithm. 
The  first  uses  only  selected  channels  which  are  determined 
to  be  free  of  any  cloud  contamination.  Then,  a  second  pass 
is  run  using  all  available  channels  not  considered  to  be 
seriously  contaminated  by  clouds.  The  determination  of 
channel  contamination  is  discussed  in  detail  by  Smith  et 
al_.  (1986). 

There  are  several  differences  between  the  "old"  and 
"new"  versions  of  the  physical  algorithm.  The  "new" 
algorithm  solves  for  temperature  and  water  vapor 
concurrently ,  whereas  the  two  parameters  are  calculated 
separately  of  each  other  in  the  "old"  scheme. 
Consequently,  any  error  in  one  parameter  would  be  carried 
over  to  the  other  in  the  Smith  (1983)  procedure.  Smith  et_ 
aJ -  (1986)  noted  that  the  simultaneous  solution 
"  alleviates  the  problem  of  the  interdependence  of  the 
radiation  observation  upon  the  three  parameters" 
(temperature,  water  vapor,  and  surface  skin  temperature) . 
Another  significant  difference  concerns  the  initial  water 
vapor  profiles  used  by  the  algorithms.  The  "new"  version 


utilizes  a  first  guess  profile  rather  than  a  50V.  RH 
assumption  based  on  the  initial  temperature  sounding. 
Since  a  first  guess  water  vapor  profile  should  be  superior 
to  a  crude  502  RH  approximation,  improved  water  vapor 
values  would  be  expected  from  the  “new"  algorithm. 

Documentation  for  the  simultaneous  algorithm  is 
limited;  however.  Dr.  C.  M.  Hayden  of  NESDIS  (personal 
communication,  1987)  noted  that  it  should  produce  little 
change  in  temperature  performance.  On  the  other  hand, 
considerable  improvement  is  expected  in  humidity. 

Even  with  advancements  in  the  retrieval  process  over 
the  past  few  years,  there  are  still  many  possible  sources 
of  error  in  the  "flow  chart"  from  the  beginning  step  of 
measuring  radiation  from  the  space  platform  itself  to  the 
final  output  of  temperature  and  dewpoint.  For  example,  one 
must  consider  the  biases  and  noise  of  the  radiometer. 
Also,  the  VAS  instrument  is  a  volume  sensor,  not  providing 
a  point  source  measure  as  done  by  radiosondes. 
Specifically,  VAS  has  poor  vertical  resolution  because  of 
the  inability  to  use  spectral  wavelengths  which  have  high 
vertical  resolution  weighting  functions  (Fig.  1). 
Similarly,  radiance  measurements  are  representat i ve  of  a 
horizontal  area  as  well.  Finally,  one  must  consider  the 
retrieval  algorithm  itself.  Truncation  errors  introduced 
by  the  mathematics  of  matrix  inversion  and  the 


introduction  of  a 


"first  guess”  profile  are  additional 


concerns  o-f  which  the  -forecaster  must  be  aware.  With  the 
above  mentioned  limitations  in  mind,  an  examination  of  a 


case  in  July  1982  is  presented. 


Since  one  of  the  goals  was  to  investigate  changes  in 
performance  between  results  of  the  "new"  and  "old" 
physical  algorithms,  the  case  studied  by  Fuelberg  and  Funk 
<1987)  was  chosen.  Consequently,  sounding  times  and  areal 
coverage  were  similar  to  those  of  their  research.  VAS 
soundings  were  available  over  the  central  United  States  at 
5  times:  1300,  1600,  1700,  2000,  and  2300  GMT  21  July 
1982.  The  satellite  operated  in  the  DI  mode  for  the  1300 
and  1600  GMT  measurements,  but  in  the  DS  mode  for  the 
latter  3  times. 

VAS  retrievals  were  prepared  at  the  University  of 
Wisconsin  on  the  Man-computer  Interactive  Data  Access 
System  <  Me  I  DAS)  (Smith  et_  al_.  ,  1979)  .  An  NMC  12  h  forecast 
valid  1200  GMT  was  used  as  the  "first  guess”  profile 
for  the  1300  GMT  retrieval  run.  Thereafter,  previous  VAS 
profiles  were  employed.  Retrieval  data  were  available  at 
950  mb  as  well  as  the  "mandatory"  levels  of  the  radiosonde 
code,  i.e. ,  the  surface,  and  850,  700,  500,  300,  250,  200, 
150,  and  100  mb. 

McIDAS  was  used  to  edit  any  profiles  contaminated  by 
clouds  or  inconsistent  with  surrounding  data.  Then,  after 


check  was  performed.  Several  deletions  were  done 
subjectively  by  the  author,  based  on  obvious  data 
inconsistencies.  Overall,  the  amount  o-f  editing  was  small, 
as  only  a  few  retrievals  were  eliminated.  The  number  of 
soundings  available  after  editing  is  given  in  Table  2  for 
each  time.  The  midwestern  United  States  was  the  area  of 
interest,  and  data  coverage  there  was  excellent  with  only 
a  few  data  gaps  due  to  cloud  cover  being  evident  (Fig.  3). 
The  average  spacing  for  the  retrievals  was  125  km. 

Radiosonde  data  for  1200  GMT  21  July  and  0000  GMT  22 
July  1982  were  obtained  from  the  National  Climatic  Center. 
Copies  of  NMC  facsimile  charts  were  gathered  from  local 


Table  2.  Number  of  soundings  within  the  computational 
domain  after  editing. 


Sounding  Time 
(GMT) 


Number  of 
soundings 


1300  119 

1600  127 


1700 


113 
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5.  Synoptic  conditions 

Conditions  on  21  July  1982  began  much  like  those  of 
the  previous  few  days.  The  main  baroclinic  zone  stretched 
along  the  northern  border  of  the  United  States,  while 
upper  level  dynamics  elsewhere  were  weak. 

The  1200  GMT  NMC  surface  chart  <Fig.  4)  depicted  a 
weak  stationary  front  extending  from  a  low  pressure  area 
over  eastern  Colorado  through  central  Missouri  and 
northern  Georgia.  Another  cold  front  reached  from 
northeastern  Minnesota  to  northern  Nebraska  and  central 
Montana,  with  a  slightly  cooler  and  drier  airmass  from 
Canada  to  its  north.  Skies  over  several  areas  contained 
remnants  of  nighttime  convection  (not  shown).  Showers  and 
thunderstorms  had  persisted  through  the  night  over  most  of 
Minnesota,  the  eastern  half  of  Nebraska,  and  northern 
Missouri.  In  addition,  scattered  thunderstorms  were  in 
progress  over  most  of  the  Gulf  Coast.  Fog  was  evident  over 
most  of  West  Virginia,  eastern  Kentucky,  and  northern 
Louisiana.  No  apparent  cloud  features  were  associated  with 
the  front  from  Nebraska  into  Montana.  Infrared  imagery 
(not  shown)  indicated  high  level  clouds  covering  South 
Carolina  and  Georgia,  most  likely  debris  from  convection 


of  the  previous  night. 
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4.  Analyses  at  1200  GMT  21  July  1982.  At  the 
sur-face,  isobars  are  labeled  in  the  usual 
convention.  At  350  mb,  heights  are  solid 
where  45  represents  1450  m,  and  isodroso— 
therms  (dashed)  are  in  "C.  At  500  mb, 
heights  are  solid  where  591  denote  5910  m 
and  temperatures  (dashed)  are  in  *C. 


At  850  mb  (Fig.  4> ,  winds  generally  were  light  and 
variable  over  most  of  the  southeastern  United  States. 
Wind  plots  (not  shown)  revealed  convergence  along  the 
front  over  the  Northern  Plains,  although  there  was  little 
convergence  along  the  Ohio  Valley  boundary.  A  moisture 
axis  (dewpoints  >  15  *C)  extended  from  North  Dakota 
through  the  middle  Mississippi  River  Valley. 

RAOB— derived  Lifted  Indices  (LI)  at  1200  GMT  are  given 
in  Fig.  5.  Most  of  Kansas,  Missouri,  Arkansas,  western 
Tennessee,  and  northern  Mississipi  were  relatively 
unstable,  with  values  less  than  2.  A  smaller  instability 
center  was  located  over  central  Minnesota.  On  the  other 
hand,  the  Great  Lakes  region  was  stable,  with  values 
greater  than  10. 

At  500  mb  (Fig.  4>  a  strong  developing  trough 
stretched  from  southern  Minnesota  through  central 
Missouri.  Twelve  hour  height  falls  of  40  m  could  be  found 
to  its  east.  Strangely,  however,  the  thermal  analysis  did 
not  indicate  a  pronounced  cold  tongue  behind  or  coincident 
with  the  height  feature.  Ridging  extended  ahead  of  the 
trough  from  Michigan  to  Tennessee,  then  west  to  a  major 
center  over  Colorado.  Using  VAS  profiles  from  the  Smith 
(1983)  algorithm,  Fuelberg  and  Funk  (1987)  calculated 
large  scale  vertical  motions  for  this  case  and  found  an 
area  of  upward  motion  to  the  east  of  the  500  mb  trough. 

Even  higher,  at  200  mb  (not  shown) ,  a  weak  trough 


stacked  appropri atel y  with  its  500  mb  counterpart; 
however,  there  was  no  well  defined  speed  maximum  nearby. 
With  the  exception  of  the  weak  short  wave,  flow  generally 
was  zonal  over  the  entire  region.  Winds  were  stronger 
upwind  of  the  trough,  suggesting  that  further 
amplification  was  likely. 

By  mid-morning,  i.e.,  1600  GMT  (Fig.  6),  the  cirrus 
debris  over  the  Southeast  had  dissipated.  Animated 
sequences  of  visible  images  revealed  an  apparent 
circulation  over  the  Northern  Plains  which  was  moving 
southeastward  into  northeastern  Iowa.  Showers  associated 
with  the  circulation  drifted  slowly  to  the  east  into 
Wisconsin.  Elsewhere,  showers  and  thunderstorms  had 
developed  over  the  northern  half  of  Florida. 

The  1800  GMT  NMC  surface  chart  (Fig.  7)  showed  a  low 
pressure  system  over  eastern  Iowa,  confirming  the 
satellite-detected  circulation  center  (Fig.  6).  Attached 
to  this  surface  low  was  a  warm  front  over  the  Ohio  River 
Valley  which  had  earlier  been  classified  as  stationary.  A 
surface  streamline  analysis  (not  shown)  performed  on  the 
Florida  State  University  Department  of  Meteorology’s  KLIMA 
(Kindly  Local  Interactive  Meteorological  Access)  system 
revealed  well  defined  convergence  around  the  surface  low. 
Similarly,  KLIMA-der i ved  analyses  of  surface  moisture 
convergence  indicated  strong  influx  in  the  vicinity  of  the 
surface  low  and  southeastward  along  the  warm  front. 


Fig.  6.  Hourly  visible  sate 
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During  the  afternoon  hours,  visible  satellite  imagery 
(Fig.  6)  showed  typical  "popcorn"  cumulus  developing  over 
most  of  the  Southeast.  Strong  convection  continued  over 
northern  Florida  and  southern  Georgia,  while  clusters  o-f 
new  thunderstorms  were  developing  over  western  Kentucky, 
central  Arkansas,  and  southern  Illinois  by  1800  GMT.  The 
latter  two  regions  appeared  to  -form  along  outflow 
boundaries  from  the  nighttime  thunderstorms  in  Missouri.  A 
strong  line  of  thunderstorms  developed  over  southwestern 
Illinois  (apparently  along  the  warm  front  noted  at  1800 
GMT  (Fig.  7)).  The  NSSFC  issued  a  severe  thunderstorm 
watch  over  most  of  Illinois  for  this  developing  line,  and 
radar  summaries  (not  shown)  indicated  tops  exceeding  18  km 
during  the  late  afternoon.  Visible  and  infrared  imagery 
suggested  that  the  500  mb  trough  had  formed  a  cut-off  low, 
with  several  small  vorticity  lobes  rotating  around  it. 

By  0000  GMT  22  July  1982,  the  NMC  surface  analysis 
(Fig.  8)  indicated  that  the  surface  low  had  moved  to 
northern  Illinois.  The  cold  front  over  the  Plains  had 
pushed  southward  to  central  Missouri,  while  a  squall  line 
was  depicted  over  eastern  Illinois.  At  850  mb,  the 
afternoon  convection  had  distorted  parts  of  the  moisture 
field;  however,  the  squall  line  continued  to  move  with  the 
moisture  axis.  Winds  still  were  light  and  variable  over 
most  of  the  region,  with  only  weak  convergence  noted  along 


the  Plains  boundary. 


The  0000  GMT  LI  analysis  (Fig.  5)  revealed  that  an 
area  of  instability  had  developed  over  central  Nebraska 
and  Kansas.  Meanwhile,  two  other  smaller  instability 
centers  were  located  over  central  Wisconsin  and  northern 
Alabama.  The  Great  Lakes  remained  relatively  stable,  as 
values  continued  to  exceed  10. 

The  500  mb  analysis  (Fig.  8)  indicated  that  the  trough 
had  indeed  become  cut-off,  with  the  center  just  southwest 
of  Green  Bay,  Wisconsin,  and  the  trough  trailing  southward 
into  central  Illinois.  The  thermal  field  now  was  better 
defined,  with  a  cold  -10*C  pocket  located  in  the  vicinity 


of  the  low.  Conditions 

at 

200 

mb 

(not 

shown ) 

consistent  with  those  at 

500 

mb 

and 

those 

seen 

earlier  (Fig.  4). 

In  summary,  21  'July  was  a  "typical"  summer  day  for 
severe  weather,  with  reports  of  one  tornado,  two  cases  of 
large  hail,  and  several  instances  of  damaging  wind 
reaching  the  NSSFC  by  late  evening.  Most  of  the  severe 


activity  was  confined  to  the  state  of  Illinois. 


One  must  understand  the  attributes  and  limitations  o-f 
VAS  soundings  before  using  them  to  calculate  derived 
atmospheric  parameters.  For  example,  do  the  data  exhibit 
spatial  and  temporal  continuity,  i.e. ,  is  there  good 
relative  agreement  between  successive  times  and  locations? 
Also,  what  are  the  magnitudes  of  bias  errors,  and  how 
great  is  the  "noise"  (RMS  errors)?  Three  approaches  will 
be  used  to  answer  these  questions.  First,  representati ve 
data  plots  will  be  examined,  followed  by  point 
comparisons  between  VAS-derived  and  RAOB-derived 
soundings.  Finally,  using  procedures  similar  to  those  of 
Fuelberg  and  Funk  (1987),  area  averaged  temporal 
characteri sti cs  will  be  presented. 

For  the  purpose  of  retrieval  assessment,  temperature, 
dewpoint,  and  height  data  provided  by  the  national  RAOB 
network  are  assumed  "perfect",  i.e.,  they  are  considered 
to  be  the  true  atmospheric  values  at  the  time  and  location 
of  the  sensor.  Of  course,  the  radiosonde  is  not  a 
"perfect"  measuring  device  (Lenhard,  1970) ,  and  some  of 
the  VAS  errors  assumed  from  the  calculations  that  follow 
actually  will  be  attributable  to  this  imperfect  standard. 
Furthermore,  as  stated  earlier,  the  sonde  is  a  "point" 


measuring  device  while  the  VAS  instrument  is  a  volume 


sensor . 

Figs.  9  and  10  depict  the  evolution  o-f  VAS-derived  850 
mb  dewpoint  and  500  mb  temperature,  respectively.  The 
author  performed  the  subjective  analyses  using  standard 
meteorological  techniques  learned  as  a  severe  weather 
forecaster  at  the  Air  Force  Global  Weather  Central 
(AFGWC) .  Where  possible,  ridges,  troughs,  moisture  axes, 
etc.,  were  enhanced  without  conflicting  with  the  data, 
i.e. ,  al 1  data  were  retained  in  the  analysis.  As  can  be 
seen,  both  isotherms  and  i sodrosotherms  are  relatively 
smooth  with  no  discontinuities. 

The  VAS  humidity  patterns  at  850  mb  (Fig.  9)  contain 
several  interesting  features.  The  moisture  axis  initially 
stretching  from  southern  Minnesota  to  the  central  Gulf 
Coast  has  reasonable  continuity.  Specifically,  there  is  a 
slight  eastward  progression  throughout  the  afternoon 
hours,  and  its  relative  strength  remains  fairly  constant 
at  16-1S°C.  Much  drier  conditions  are  found  over  the 
northeast  corner,  i.e.,  on  the  cold  side  of  the  frontal 
system  located  over  the  Ohio  River  Valley  (Figs.  4,  7, 
8).  Especially  interesting  is  the  development  of  a 
relatively  small  scale  "dry"  pocket  in  the 
Mi ssour i -I 1 1 inoi s  area  between  1600-1700  GMT.  This  feature 
will  be  examined  in  greater  detail  in  a  later  chapter. 

In  general,  the  broad  scale  aspects  of  VAS-derived 


dewpoint  patterns  at  1300  and  2300  GMT  <Fig.  9)  are 
closely  related  to  sonde-derived  versions  at  1200  and  0000 
GMT  (Figs.  4  and  8,  respectively).  Agreement  between  the 
moisture  axis  over  the  Mississippi  River  Valley  is  the 
prime  example.  It  is  clear  that  the  closer  spacing  of  the 
VAS  retrievals  has  provided  considerably  smaller  scale 


resol uti on 

of 

features  than 

is  possible 

from 

the 

radiosonde 

network.  Dn  the 

other 

hand , 

val ues 

at 

i ndi vi dual 

RA0B 

sites  sometimes 

differ 

between 

the 

two 

data  types.  This  aspect  is  explored  in  later  paragraphs. 

Although  the  satellite's  500  mb  temperature  -fields 
(Fig.  10)  are  not  as  well  defined  as  the  just  described 
850  mb  dewpoints,  several  features  do  exhibit  good 
continuity.  For  example,  a  cold  trough  initially  over  Iowa 
moves  eastward  and  out  of  the  data  domain  during  the 
afternoon.  Unfortunately,  since  it  is  located  near  a  cloud 
induced  data  gap  (Fig.  3) ,  one  cannot  have  great 
confidence  in  its  exact  location  or  minimum  temperature.  A 
second  cold  axis  consistently  located  over  Nebraska  is 
more  pronounced  by  the  ending  time.  Finally,  a  cold  pocket 
over  southern  Georgia  is  quasi -stationary.  Once  again, 
broad  scale  features  of  the  VAS  analyses  have  similarities 
to  those  from  the  standard  source  (Figs.  4,  8).  Also,  VAS 
temperatures  vary  several  degrees  from  one  time  to 
another . 

In  order  to  estimate  the  accuracy  of  VAS  measurements, 


retrievals  were  paired  with  nearby  r adi osonde-der i ved 


"ground  truth"  data.  Then,  the  VAS-derived  values  were 
subtracted  -from  the  sonde-derived  measurements.  This 
procedure  is  similar  to  that  employed  by  previous  authors, 
e.g.,  Phillips  aj (1979),  Gruber  et  al _.  (1982).  Three 
criteria  had  to  be  met  -for  each  comparison.  First,  the  VAS 
pro-file  being  paired  had  to  be  within  a  100  km  radius  o-f 
the  RAOB  location.  Second,  the  100  km  criterion  had  to  be 
met  at  both  the  1200  and  0000  GMT  RAOB  sounding  times. 
Third,  the  same  RAOB  sites  had  to  be  available  in 
evaluating  both  the  simultaneous  and  the  physical 
retrievals,  i.e.,  the  same  numbers  and  locations  of  RAOBs. 
There  are  at  least  two  unavoidable  limitations  to  the 
pairing  procedure.  First,  although  the  100  km  value  was 
required  to  obtain  a  sufficient  number  of  pairings, 
structure  function  analyses  of  radiosonde  and  VAS  data 
have  showed  that  appreciable  horizontal  gradients  still 
exist  at  scales  smaller  than  100  km,  especially  for  water 
vapor  (Fuelberg  and  Meyer,  1985).  Second,  there  was  a  0-2 
h  difference  between  the  two  data  sources,  i.e.,  1200  GMT 
RAOBs  (released  near  1100  GMT)  were  compared  to  1300  GMT 
VAS  profiles,  and  0000  GMT  RAOBs  (released  near  2300  GMT) 
were  compared  to  2300  GMT  VAS  locations.  Thus,  this  lack 
of  complete  spatial  and  temporal  agreement  will  be  a 
source  for  some  of  the  calculated  discrepancies  that  are 


described  next. 


Table  3  gives  results  for  both  the  simultaneous  and 
older  physical  retrievals  (used  originally  by  Fuelberg  and 


Funk,  1987).  Seven  RAOB-OAS  data  pairs  met  the  three 
criteria.  These  relatively  small  sample  sizes  are  a 
limitation  to  the  evaluation.  It  was  possible  to  obtain  16 
RAOB-VAS  pairings  for  the  simultaneous  retrievals  by 
dropping  the  third  criterion;  however,  those  results  (not 
shown)  are  similar  to  those  of  the  smaller  sample  size 
(Table  3).  For  comparison  purposes,  previously  published 
results  for  the  older  algorithm  (Jedlovec,  1985)  and  the 
new  algorithm  (Wade  et_  ai. .  ,  1985)  also  are  included  in  the 
Table. 

i 

I  The  simultaneous  (new)  algorithm  appears  to  provide 

better  estimates  of  water  vapor  profiles  (Table  3); 
however,  the  temperature  soundings  are  inferior  to  those 
of  the  older  approach.  Specifically,  mean  dewpoint  bias 
errors  for  the  simultaneous  retrievals  are  considerably 
smaller  than  those  from  the  physical  algorithm.  Both 

I  procedures  yield  dewpoints  that  are  warmer  (positive  bias) 

'  than  RAOBs  in  the  lower  troposphere,  but  colder  aloft.  RMS 

errors  for  the  simultaneous  retrievals  are  comparable  to 

i 

I  those  of  the  older  approach  at  850  and  700  mb,  but  the 

newer  results  clearly  are  superior  at  the  higher  levels. 
Current  simultaneous  mean  bias  errors  are  smaller  than 

I  those  of  Jedlovec  (1985),  but  there  is  no  consistent 

relation  with  his  RMS  errors.  In  addition,  simultaneous 

I 


Table  3.  A  comparison  of  OAS  retrievals  paired  with  near — 
by  RAOBs.  The  le-ftmost  value  o-f  each  entry  is  the 
mean  error  while  the  rightmost  value  is  RMS 
error.  Units  are  *C  -for  dewpoint  and  temperature 
and  m  -for  height. 


Si mul taneous 

Physical 

Pressure 
level  (mb) 

This 

study 

Wade 

et  al  . 

( 1985) 

This 

study 

Jedl ovec 
( 1985) 

300 

400 

-3. 6/4. 4 
-2. 2/4. 4 

Dewpoint 

-6. 5/8. 3 
-7. 6/6. 3 

4. 0/3. 6 

500 

0.8/4. 4 

-9. 2/5. 4 

-2. 4/4.6 

700 

2. 0/4. 7 

4. 9/4. 3 

-5. 9/5. 5 

S50 

0.3/3. 3 

0.4/7. 8 

4.3/3. 1 

-1.8/2. 3 

200 

-4. 0/1. 5 

Temperature 

-1.9/1. 7 

-0. 1/1.3 

300 

-3. 3/1. 7 

-1. 5/2.0 

-1.7/1. 6 

400 

-2. 9/1. 5 

-2. 0/2. 3 

-1.7/1. 4 

500 

-1.5/1. 7 

0.3/1. 8 

-1.5/1. 3 

-1.7/1. 6 

700 

1.9/2. 2 

0.2/2. 2 

0.4/1. 1 

850 

2. 2/2.0 

0.9/2. 3 

0. 5/2.0 

-2. 5/1. 5 

200 

-42. 8/37.5 

Hei ght 

-17. 1/37.3 

-42.0/21.0 

300 

-9.4/36.3 

1.6/23.4 

-31.0/21.0 

400 

9.3/36.7 

3.  1/16.9 

-13.0/18.0 

500 

19.2/34.8 

10. 1/16.2 

-4.8/12.0 

700 

9.4/18.9 

8.9/  9.7 

4.2/  7.9 

850 

-3.0/10.3 

3.4/  5.5 

4.3/  6.8 
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errors  at  850  mb  are  smaller  than  those  of  Made  et  al 


(1985)  . 


The  greatest  dewpoint  difference  (simultaneous 


retrievals)  at  a  single  level  occurred  over  Denver  at  300 


mb  (-20.6*0.  Fig.  11  (12)  displays  the  "best"  ("worst") 


dewpoint  profile,  i.e.,  considering  the  entire  sounding. 


Errors  averaged  2.0  *C  at  each  level  for  the  "best" 


sounding,  but  exceeded  10*C  between  the  surface  and  400  mb 


for  the  "worst"  case.  The  reduced  vertical  resolution  of 


the  VAS  retrievals  is  clearly  evident.  Furthermore,  it  is 


interesting  that  VAS-derived  dewpoint  retrievals  having 


good  agreements  with  the  RAQB  versions  often  correspond  to 


temperature  profiles  differing  greatly  from  those  of  the 


sondes,  and  vice  versa.  As  an  extreme  example,  the  pairing 


which  exhibits  the  best  overall  dewpoint  agreement  (Fig. 


11)  also  has  the  worst  temperature  comparison. 


Mhen  temperatures  from  the  two  algorithms  are  compared 


with  RAOBs  (Table  3),  the  simultaneous  version  exhibits 


increased  magnitudes  of  mean  error.  The  new  soundings  tend 


to  be  slightly  warmer  than  RAOBs  in  the  lower  troposphere: 


however,  the  bias  becomes  negative  (cooler)  by  the  500  mb 


level,  increasing  to  -4.0*C  at  200  mb.  This  tendency  is 


consistent  with  that  of  the  older  physical  retrievals.  RMS 


errors  vary  little  with  altitude  and  are  approximately  1.7 


C  for  both  sets  of  soundings.  The  greatest  temperature 


error  at  a  single  level  also  occured  in  the  Denver  profile 
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at  300  mb  (+8.4  *C).  Fig.  13  (11)  shows  the  "best" 
("worst")  VAS/RAOB  temperature  comparisons.  Errors  at  each 
level  are  1.0*C  (  4.0*0  from  400  mb  to  the  surface  for 


the  "best"  ("worst")  sounding. 

Since  heights  are  obtained  by  integrating 
temperature  and  dewpoint  through  the  hypsometric  equation, 
biases  for  height  correspond  to  those  of  the  temperatures 
(Table  3).  Greatest  positive  bias  occurs  at  500  mb  in 
both  cases,  but  values  rapidly  become  negative  at  higher 
levels.  Greatest  negative  biases  occur  at  200  mb.  The 
increasingly  large  negative  values  in  the  upper 
troposphere  result  from  the  cold  bias.  RMS  errors  for  the 
simultaneous  retrievals  often  are  twice  those  from  the 
physical  version. 

Although  the  newer  algorithm  produces  smaller  mean  and 
RMS  dewpoint  errors,  it  is  disconcerting  that  the  mean 
temperature  errors  are  greater.  Dr.  C.  M.  Hayden  of  NESDIS 
has  kindly  provided  an  explanation  for  this  phenomenon 
(personal  communication,  1987).  Specifically,  the  VAS 
radiometer  now  is  known  to  produce  biased  radiances, 
possibly  due  to  temperature  variations  on  its  filter 
wheel.  To  control  the  mean  retrieval  errors  that  result 
from  this  input,  NESDIS  alters  the  radiances.  These 
radiance  corrections  had  been  updated  weekly,  but  future 
plans  call  for  daily  updates.  Based  on  the  increased  mean 
temperature  bias  (Table  3),  it  is  assumed  that  radiance 


corrections  applied  when  the  simultaneous  retrievals  Mere 
prepared  in  1985  are  not  the  same  as  those  used  in 
creating  the  older  physical  version  in  1982.  This  aspect 
of  the  retrieval  process  underscores  the  importance  of 
emphasizing  VAS-derived  patterns  at  single  times  and 
pattern  evolutions  over  short  periods.  Conversely,  one 
should  not  stress  magnitudes  of  VAS  products  at  individual 
locations  or  their  comparison  with  RAOBs.  This  aspect  has 
been  noted  operationally  by  personnel  at  NSSFC  (Anthony 
and  Wade,  1983) . 

Fuelberg  and  Funk  (1987)  calculated  mean  soundings 
Mi  thin  several  geographic  "boxes"  over  the  midMestern 
United  States  in  order  to  describe  temporal  variations  in 
the  physical  retrievals.  To  compare  performance  changes 
between  the  two  algorithms,  similar  boxes  were  chosen 
(Fig.  14)  for  this  study.  Box  "A"  enclosed  a  large  area 
which  was  relatively  cloud-free  during  the  entire  day, 
while  Box  "B"  was  centered  over  an  area  which  experienced 
strong  afternoon  convection.  Using  simple  averaging,  a 
mean  profile  was  created  for  each  box  at  each  time.  For 
Box  A,  approximately  30  soundings  from  the  new  algorithm 
were  averaged  at  each  time,  whereas  there  were  40  from  the 
older  algorithm.  The  averages  for  Box  B  consisted  of 
approximately  15  (20)  soundings  for  each  time  from  the  new 
(old)  data.  To  clearly  reveal  the  temporal  tendencies, 


area-averaged  profiles 


at 


consecutive 


t  i  mes 


were 


only 


subtracted  -from  each  other.  For  the  sake  of  brevity, 
results  -from  the  1300-1600  GMT  (Fig.  15)  and  2000-2300  GMT 
(Fig.  16)  calculations  are  presented  here.  For  comparison 
purposes,  results  -from  Fuelberg  and  Funk  (1987)  also  are 
plotted. 

Fuelberg  and  Funk  (1987)  noted  that  physically-derived 
temperatures  -from  near  the  surface  to  approximately  400  mb 
exhibited  strong  warming  during  the  morning  (Fig.  15),  but 
there  was  compensating  cooling  at  higher  levels.  These 
trends  reversed  during  the  afternoon  (Fig.  16) .  The 
variations  could  not  be  explained  from  temperature 
advection,  and  one  would  not  suspect  daily  surface  heating 
to  affect  levels  much  above  the  boundary  layer. 
Therefore,  they  mostly  attributed  the  unrealistic  trends 
to  deficiencies  in  the  retrieval  algorithm.  With  the  newer 
soundings,  the  simultaneous  algorithm  performs 
considerable  better  during  the  afternoon  (Fig.  16); 
however,  there  still  are  major  variations  during  the 
morning  (Fig.  15)  that  appear  to  be  algorithm  related. 

The  area  averaged  dewpoint  profiles  of  Fuelberg  and 
Funk  (1987)  exhibited  tendencies  similar  to  those  of 
temperature  (Figs.  15,  16).  That  is,  morning  values 

increased  (decreased)  in  the  lower  (upper)  troposphere, 
and  there  was  a  reversal  between  morning  and  afternoon. 
Since  there  was  no  meteorological  explanation  for  the 
variations,  it  appeared  that  humidity  calculations  were 
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closely  related  to  the  questionable  temperatures.  With  the 
newer  simultaneous  retrievals  (Figs.  15,  16),  unrealistic 
trends  remain.  Specifically,  there  is  a  large  increase  (3 
—  4  *C)  in  mid— level  moisture  during  the  morning  hours 
followed  by  drying  in  the  afternoon,  but  again  this  cannot 
be  substantiated  by  satellite  imagery  or  conventional 
data. 

In  summary,  from  the  standpoint  of  accuracy,  the  new 
simultaneous  algorithm  should  be  rated  better  than  the 
Smith  (1983)  scheme.  Moisture  values  are  closer  to  the 
"truth",  although  temperature  accuracy  is  slightly  worse. 
Diurnal  variations  have  been  reduced,  but  not  eliminated, 
and  some  of  the  compensating  effects  have  been  removed. 
Armed  with  this  information  about  the  basic  data,  several 
different  thermodynamic  parameters  related  to  severe  storm 
forecasting  will  be  presented  in  the  following  chapters. 
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7.  Non-i ntegr ated  stability  parameters 

Summer  is  the  most  difficult  season  for  forecasting 
the  location  and  intensity  of  convection  because  the 
atmosphere  generally  is  sufficiently  unstable  for  even 
minor  ascent  ("trigger"  mechanisms)  to  initiate  possible 
severe  storm  activity.  A  meteorologist  needs  to  know  the 
patterns  of  atmospheric  stability  together  with  a 
knowledge  of  the  vertical  motion.  Synoptic  scale  and 
mesoscale  "triggers"  on  21  July  were  described  in  Chapter 
5;  this  chapter  considers  stability. 

Several  thermodynamic  indexes  are  used  by  operational 
meteorologists  to  evaluate  stability.  They  usually  are 
derived  from  input  at  various  atmospheric  levels  and  can 
be  placed  into  two  distinct  groups,  integrated  and 
non-i ntegrated .  The  integrated  variety  is  based  on 
numerous  data  over  a  thick  layer,  rather  than  measurements 
at  a  few  discrete  levels  (the  non-i ntegrated  type). 
Examples  of  integrated  parameters  include  precipitable 
water  and  buoyant  energy.  This  chapter  will  deal  with  the 
non-i ntegr ated  parameters,  specifically,  the  Lifted  Index 
(LI),  K  Index  (KI),  and  Total  Totals  (TT)  Index. 


a.  Lifted  index 
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The  LI  is  the  most  widely  used  stability  index  of  the 
National  Weather  Service  and  is  calculated  every  12  h  from 
all  RAOBs  in  the  United  States.  The  LI  is  based  on  a 
temperature  comparison  between  a  saturated  parcel  which 
has  been  lifted  adi abati cal 1 y  to  500  mb  and  the  actual  500 
mb  sounding  temperature.  The  parcel  temperature  is 
subtracted  from  the  sounding  temperature  to  produce  the 
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index  value.  When  the  parcel  is  warmer  than  the 
environment  (negative  index),  the  atmosphere  is  considered 
unstable.  An  atmosphere  with  an  LI  of  -4  or  less  is 
classified  as  having  severe  storm  potential.  Temporal 
variations  in  LI  also  provide  useful  information.  Thus, 
the  NSSFC  calculates  hourly  changes  in  LI  by  blending 
surface  data  with  the  12  h  RAOBs  (Hales  and  Doswell, 
1982) ,  thereby  helping  forecasters  pinpoint  areas  of 
severe  weather  potential. 

There  are  a  variety  of  schemes  for  calculating  Lis, 
with  differences  mainly  due  to  the  method  of  initially 
specifying  the  parcel.  The  LI  values  presented  in  this 
paper  ware  based  on  the  following  algorithm:  First,  the 
mean  temperature  and  dewpoint  were  calculated  for  the 
lowest  kilometer  by  using  simple  averaging  of  all 
retrieved  values  in  that  layer.  One  may  argue  that  this 
negates  inclusion  into  the  non-integration  category; 
however,  the  layer  is  much  more  shallow  than  employed  in 
tne  integrated  variety  described  later.  After  mean  values 
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were  determined,  the  Lifting  Condensation  Level  (LCL)  was 
computed  by  raising  the  parcel  dry  adi abatical 1 y  until 
saturation.  Then,  the  parcel  was  cooled  moist 
adi abati cal 1 y  to  500  mb  where  the  temperature  comparison 
was  made. 

Fig.  17  displays  fields  of  VAS-derived  LI  for  the  five 
observation  times.  As  with  the  dewpoint  analyses  described 
earlier  (Fig.  9),  the  Lis  have  good  spatial  and  temporal 
continuity.  Several  features  in  the  1300  GMT  field  undergo 
interesting  evolutions  that  should  be  discussed  in  detail. 
First,  the  relative  minimum  (unstable)  in  eastern  Missouri 
(-6)  slightly  loses  intensity  during  the  day,  becoming  as 
stable  as  -3.  The  area  becomes  a  well  defined  relative 
maximum  (more  stable)  by  2000  GMT  and  continuing  to  2300 
GMT  mostly  because  the  surrounding  areas  become  much  less 
stable.  The  comparitive  afternoon  stability  over  eastern 
Missouri  appears  due  to  early  morning  convection  in  the 
area.  Specifically,  a  strong  outflow  boundary  was  evident 
in  satellite  imagery  (Fig.  6),  suggesting  drier,  downdraft 
air  from  the  early  morning  storms.  This  is  reflected  in 
the  850  mb  dewpoint  fields  at  1700  GMT  (Fig.  9)  by  a 
relative  minimum  (-10*0  over  eastern  Missouri.  This  drier 
region  at  850  mb  is  consistent  with  that  from  routine 
surface  data  (not  shown).  Thus,  VAS  suggests  that  the 
apparent  downdraft  air  is  not  confined  merely  to  the 


boundary  layer. 
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To  further  explore  the  apparent  outflow,  VAS 
retrievals  for  the  area  are  given  in  Fig.  IS;  actual 
sounding  locations  are  noted  in  Fig.  17  by  the  symbol  0. 
Profiles  east  and  west  of  the  boundary  at  2000  GMT  (Fig. 
18a)  reveal  marked  drying  from  the  ground  to  the  middle 
troposphere  west  of  the  boundary.  Examination  of  profiles 
at  1300  and  2300  GMT  (Fig.  18b)  reveals  a  large  daytime 
decrease  in  dewpoint  from  the  surface  to  near  6>0u  mb.  The 
slight  increase  in  moisture  in  the  middle  troposphere 
might  be  a  result  of  the  convection,  i.e. ,  debris.  The  dry 
low  level  air  would  limit  further  storm  development  over 
the  area,  and  this  is  confirmed  by  satellite  imagery 
indicating  no  convection  over  Missouri  after  1700  GMT. 
Strong  diurnal  warming  is  noted  near  the  surface,  however 
the  2-3  *C  increases  above  850  mb  are  suspect. 

Another  evolving  feature  at  1300  GMT  (Fig.  17)  is  the 
relative  minimum  (-6)  in  northern  Mississippi  which  seems 
to  migrate  northward  throughout  the  day.  Values  reach  -8 
at  1700  GMT,  and  -10  by  2000  GMT,  before  diminishing  at 
2300  GMT  due  to  the  storms  over  Illinois.  As  noted  earlier 
(Chapter  5)  this  was  the  area  of  the  most  intense 
afternoon  convection  (Fig.  6).  To  explore  the  rapidly 
changing  stability,  VAS-derived  soundings  for  1300  and 
2000  GMT  over  southern  Illinois  (Fig.  19)  are  examined. 
Once  again,  their  locations  are  given  in  Fig.  17,  this 


time  by  Xs.  The  Lis  for  these  soundings  change  from  —1  at 
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1300  GMT  to  -9  at  2000  GMT.  The  morning  pro-file  exhibits 
the  classic  severe  weather  pattern,  i.e.,  a  shallow  moist 
layer  near  the  surface  capped  by  an  inversion  and  drier 
air  aloft.  In  the  afternoon  sounding,  VAS  indicates  that 
surface  diabatic  heating  has  destroyed  the  morning 
inversion,  while  moisture  has  increased  from  the  surface 
to  near  400  mb.  With  the  exception  of  the  boundary  layer, 
temperatures  are  relatively  static  during  the  7  h  time 
period;  however,  slight  cooling  is  suggested  above  500  mb, 
thereby  agreeing  with  the  RA0B  analyses  (Figs.  4,  S) .  The 
1700  and  2000  GMT  LI  fields  (Fig.  17)  clearly  provide  the 
forecaster  with  excellent  guidance  since  the  strong 
(severe)  convection  forms  in  the  minimum  areas  (-8) 
defined  by  VAS  from  southern  Illinois  to  western 
Tennessee. 

The  morning  and  a-  ter noon  convection  over  southern 
Georgia  and  northern  Florida  (Fig.  6)  occurs  near  another 
LI  minimum,  -6  at  1300  GMT  (Fig.  17).  This  region  shifts 
slowly  toward  the  southwest,  while  exhibiting  increased 
instability  (-10).  Profiles  at  1300  and  2300  GMT  for 
southwestern  Georgia  (Fig.  20)  illustrate  the  changing 
conditions  (#s  in  Fig.  17).  With  these  profiles,  the  LI 
changes  from  -4  at  1300  GMT  to  -7  at  2300  GMT.  In  the 
boundary  layer,  surface  heating  warms  the  surface  (850  mb 
level)  approximately  7  *C  (3.5  °C) ,  although  dewpoints 
change  little.  During  the  morning,  VAS  indicates  that  the 


atmosphere  is  relatively  moist  -from  the  surface  to  near 


400  mb,  but  later,  with  the  exception  of  700  mb,  drying 
takes  place  at  all  levels.  The  increased  temperatures 
between  850  and  500  mb  cannot  be  explained,  and  are 
considered  suspect.  This  is  a  specific  example  of  the 
excessive  diurnal  variations  produced  by  the  simultaneous 
algorithm.  One  should  recall,  however,  that  only 
temperatures  in  the  boundary  layer  and  at  500  mb  affect 
the  LI. 

In  addition  to  the  LI  centers,  several  stability  axes 
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also  provide  useful  information.  At  1600  and  1700  GMT 
(Fig.  17),  for  example,  there  is  an  unstable  axis  (-6) 
across  northern  Arkansas  and  eastern  Oklahoma.  One  to  two 
hours  later,  satellite  imagery  (Fig.  6)  reveals  the 
development  of  a  strong  line  of  thunderstorms  in  the  area 
that  apparently  was  induced  by  an  outflow  boundary  which 
propagated  southward  from  a  morning  position  over  central 
Missouri.  As  a  second  example,  there  is  an  axis  of 
relatively  stable  air  at  all  five  times  from  southern 
Louisiana  to  southern  Mississippi.  In  this  case,  satellite 
imagery  reveal s  the  area  to  be  largely  convection  free 
during  the  entire  day.  The  1300  and  2300  GMT  VAS  soundings 
(Fig.  21)  for  this  location  (*s  in  Fig.  17)  indicate 
strong  surface  warming;  however,  dewpoints  remain  nearly 
constant.  The  relatively  small  amount  of  available 


moi sture 


apparently  is  the  reason  for  the  lack  of 
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3  retrievals  over  southern  Louisiana  at 
GO  (thin)  and  2300  GMT  (bold).  Sounding 
locations  are  given  by  the  symbol  *  ir 
rig.  17. 


convective  development-  Once  again,  the  2-3  C  increases  in 


temperature  above  the  boundary  layer  are  suspect. 

The  final  area  of  interest  is  major  instability  (-10) 
that  develops  at  2000  GMT  over  central  Nebraska  (Fig.  17). 
This  first  was  thought  to  be  a  "false  alarm",  but  radar 
charts  for  the  upcoming  evening  (not  shown)  indicate  that 
this  area  did  experience  strong  convection,  eventually 
forming  a  Mesoscale  Convective  Complex  (MCC)  over 
northeastern  Kansas.  Examination  of  1300  and  2300  GMT 
profiles  (Fig.  22)  reveals  marked  warming  near  the  surface 
(-12*0  during  the  day  and  an  increase  in  low  to  mid  level 
moisture.  (Locations  of  the  soundings  are  denoted  by  +  s 
in  Fig.  17.)  The  morning  inversion  is  exaggerated  by  VAS, 
possibly  due  to  a  problem  in  "blending"  the  surface  data 
into  the  "first  guess"  profile.  At  both  times,  however, 
the  lapse  rate  becomes  nearly  dry  adiabatic  between  350 
and  500  mb.  Although  Nebraska  is  pref erenti al 1 y  cold  at 
500  mb  (Fig.  10),  there  is  little  temperature  change  at 
these  two  sounding  sites. 

There  is  one  notable  area  where  VAS— derived  LI  (Fig. 
17)  was  not  a  good  indicator  of  convective  activity.  A 
band  of  showers  and  thunder storms  developed  over  extreme 
eastern  Texas  during  the  afternoon  (Fig.  6),  yet  LI  values 
remained  relatively  stable  (0  to  -2).  Skies  were 
cloud-free  during  the  morning  hours,  providing  abundant 


solar  heating 


and  apparently  there  was  sufficient  small 


scale  -forcing  to  "trigger"  even  this  marginally  unstable 
air. 

I 

1 

•  In  summary,  the  VAS-derived  Lis  appear  to  be  a 

,  valuable  -forecasting  tool.  That  is,  most  of  the  relatively 

unstable  areas  coinsided  with  afternoon  convection,  while 
compari t i vel y  stable  areas  were  storm  free.  Several  larger 
scale  stability  features  compared  favorably  with  their 
RAOB-derived  counterparts  (Fig.  5),  e.g.,  the  major 

|  instability  axis  over  the  Mississippi  River  Valley  and  the 

^  development  of  strong  instability  in  Nebraska.  Of  course, 

t 

VAS  was  able  to  define  these  features  at  a  finer  scale. 

i 
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I  For  example,  several  axes  of  instability  and  stability, 

|  e.g.,  the  eastern  Oklahoma  to  western  Arkansas,  and 

southern  Louisiana  locations  would  have  "fallen  through 
the  cracks"  of  the  national  RAOB  network.  One  must  notice, 
however,  that  actual  values  of  VAS-derived  LI  tend  to  be 
too  negative,  as  observed  by  Anthony  and  Wade  (1933)  and 
others.  Thus,  as  noted  earlier,  the  forecaster  should 
concentrate  on  the  rel ati ve  patterns  and  magnitudes  of  the 
features,  and  not  their  absolute  values. 

An  important  question  about  the  enhanced  resolution  of 
VAS  products  still  must  be  answered:  Do  VAS-derived  values 
provide  information  that  could  not  have  been  obtained  from 
conventional  surface  data,  i.e.,  is  VAS  merely  duplicating 
results  that  could  have  been  calculate^  from  simple 


surface  observations? 


To  answer  that  question,  Lis  (Fig. 
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23)  were  re-calculated  using  a  technique  similar  to  that 
employed  by  the  NSSFC  (Hales  and  Doswell,  1982).  The 
scheme  utilizes  500  mb  temperatures  together  with 
conventional  hourly  reports  of  temperature  and  dewpoint. 
There  was  one  difference  between  their  method  and  that 
employed  here,  i.e.,  the  current  500  mb  temperature  was 
held  constant  in  time,  whereas  NSSFC  routinely  uses  a  500 
mb  temperature  forecast  from  one  of  the  operational  NMC 
models.  Changes  in  LI  based  on  the  surface  data  variations 
then  were  calculated  by  subtracting  the  values  at  1600  and 
2000  GMT  <Fig.  23a).  A  major  assumption  of  the  Hales  and 
Doswell  (1982)  approach  is  that  the  surface  data  are 
representati ve  of  the  entire  1  km  boundary  layer.  While 
this  probably  is  valid  during  the  afternoon  hours,  it  is 
not  appropriate  during  the  morning.  Furthermore,  after 
2000  GMT,  the  surface  data  fields  were  heavily 
contaminated  by  thunderstorms ,  cooling  the  individual 
locations  by  precipitation  and  advection  processes.  Thus, 
only  results  for  the  1600-2000  GMT  time  period  are 
presented.  Correspondi ng  LI  changes  based  solely  on  VA3 
retrievals  are  given  in  Fig.  23b;  patterns  at  the 
individual  times  were  given  earlier  in  Fig.  17. 

If  VAS  merely  was  duplicating  the  surface  data,  i.e., 
not  providing  independent  information,  the  change  fields 
from  the  two  techniques  (Fig.  23)  would  be  similar. 
However,  inspection  of  the  figure  reveals  that  this  is  not 
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VAS  indicates 


the  case.  Daring  the  1600-2000  GMT  period, 
that  northwestern  Illinois,  eastern  Oklahoma,  and  eastern 
Nebraska  were  becoming  more  unstable  by  at  least  4  units. 
Some  increases  in  stability  are  noted  in  the  vicinity  of 
ongoing  convection,  e.g.,  southern  Alabama  and  central 
Georgia.  The  surf ace— derived  Lis  (Fig.  23 a)  also  display 
greater  than  -4  changes  covering  southeastern  Nebraska, 
western  Iowa,  and  western  and  southern  Missouri;  however, 
unlike  the  VAS  product,  stability  increases  occur  over 
central  Arkansas  and  eastern  Oklahoma,  and  east  of  85 °\A. 
To  further  compare  the  LI  change  fields,  Fig.  23c  shows 
differences  between  VAS-  and  surface-derived  LI  changes 
(i.e. ,  differences  between  Figs.  23a,  b) .  VAS-derived  LIs 
are  becoming  more  unstable  than  the  surf ace-deri ved  LIs 
over  Iowa  and  Illinois  as  well  as  eastern  Arkansas  (larger 
negative  values).  On  the  other  hand,  when  compared  to 
results  solely  from  the  surface  data,  VAS  indicated  that 
central  Texas,  central  Mississippi,  southern  Missouri,  and 
western  Kansas  were  becoming  either  slightly  more  stable 
or  not  as  unstable  as  surface  data  suggested. 

Based  on  Fig.  23,  it  is  obvious  that  VAS  does  not 
simply  duplicate  the  LI  variations  due  to  surface  data 
fluctuations.  Indeed,  there  were  several  large  contrasts 
in  the  patterns  of  LI  changes.  Therefore,  VAS  is  supplying 
the  forecaster  with  information  above  the  surface  that  is 
important  in  producing  stability  variations.  Using  the 


satellite  imagery  as  a  standard,  it  appears  that  the 
VAS-derived  LI  changes  are  superior  to  those  created 
solely  from  surface  data.  That  is,  areas  that  were 
becoming  more  unstable  than  indicated  by  surface  data 
often  experienced  convection  (Fig.  6),  e.g.,  eastern 
Oklahoma,  central  Illinois,  and  central  Tennessee. 
Conversely,  areas  where  VAS  indicated  more  of  an  increase 
in  stability  often  were  free  of  storms,  e.g.,  southwestern 
Missouri,  and  Louisiana,  Mississippi,  and  western  Alabama. 

A  second  question  then  arises:  What  information 
provided  by  VAS  caused  the  resulting  LI  changes  to  be 
superior?  The  two  possibilities  are  more  realistic 
information  about  500  mb  temperatures ,  i.e.,  not  the 
constant  values  assumed  in  Fig.  23,  or  improved  low-level 
moisture  depictions.  Based  on  the  sounding  evaluations 
presented  in  Chapter  £>,  along  with  analyses  of  VAS-derived 
mid-level  temperature  variability  (not  shown),  it  appears 
that  the  latter  is  the  significant  contributor. 
Specifically,  VAS ' s  ability  to  describe  moisture  fields 
both  horizontally  (Fig.  9)  and  vertically  apparently  is 
the  reason  for  its  superior  LI  patterns.  More  importantly. 
Funk  (1985)  performed  a  similar  investigation  on  Lis 
generated  from  the  physical  (Smith,  1983)  algorithm  but 
found  that  those  versions  did  not  provide  much  useful 
information  over  conventional  surface  data.  Thus,  it  is 
encouraging  that  the  new  simultaneous  profiles  are  an 


i mprovement . 

To  summarize  the  LI  findings,  their  improved 
resolution  allows  forecasters  to  concentrate  on  specific 
areas  instead  of  larger  regions.  Although  absolute  values 
often  did  not  compare  closely  with  RAOB-derived  versions, 
the  large  scale  patterns  were  similar,  and  several 
smaller — scale  features  not  apparent  in  the  RAOB  fields 
could  be  related  to  the  presence  or  absence  of  storm 
activity.  Finally,  as  just  discussed,  VAS  appears  to  give 
the  forecaster  more  than  just  a  copy  of  surface  induced  LI 
variability;  instead,  it  complements  those  data. 

b.  K  Index 

The  "K  Index"  is  a  common  parameter  for  forecasting 
summertime  airmass  convection  and  is  given  by 

K I  ~ 7  tTggo  Taoo)  "T  DPgao  (T  *700  “DP  700),  (11) 

where  T  represents  temperature,  DP  represents  dewpoint, 
and  the  subscripts  indicate  the  level  of  the  measured 
data.  The  K I  was  designed  to  estimate  thunderstorm 
potential  on  the  merits  of  thermodynamics  alone.  Thus, 
there  is  not  the  requirement  for  dynamic  lifting  that  is 
implicit  with  the  LI.  In  fact,  the  KI  does  not  perform  as 
well  as  the  LI  when  dynamic  forcing  is  present.  The  KI  is 


a  logical  choice  for  the  forecaster  of  summertime 
convection  since  large  scale  forcing  is  compar i ti vel y 
weak . 

When  radiosonde  data  are  the  input,  typical  summertime 
KI  values  range  from  20  to  50,  where  20  represents  a 
20-40/i  chance  of  a  thunderstorm,  and  50  indicates  a  1007, 
thunderstorm  probability,  along  with  a  high  likelihood  of 
a  severe  weather  event  (Miller,  1972).  Satellite  derived 
KIs  exhibit  the  same  exaggerated  values  as  the  Lis,  i.e., 
magnitudes  tend  to  be  too  large.  For  example,  KIs  for  the 
current  case  range  from  30  to  65.  No  attempt  was  made  to 
correlate  these  versions  with  the  sonde-derived  variety; 
instead,  only  relative  strengths  and  patterns  were 
examined. 

A  variation  of  the  KI  is  the  so-called  "modified"  K 
Index  (Charba,  1979)  in  which  surface  dewpoint  and 
temperature  are  averaged  together  with  the  850  mb  values. 
As  one  might  expect,  this  procedure  increases  index  values 
due  to  the  higher  averaged  temper atures.  Both  versions  of 
the  index  were  calculated  for  the  21  July  case.  Since  the 
modified  variety  consistently  outperformed  the  original 
type,  only  the  modified  KIs  (Fig.  24)  will  be  shown. 

At  1300  GMT,  the  KI  analysis  (Fig.  24)  has  maxima 
(instability)  over  eastern  Missouri  (51),  northern 
Mississippi  (50),  and  northwestern  Georgia  (50).  During 
the  morning  hours,  the  maximum  over  Missouri  appears  to 


drift  slightly  to  the  east,  generally  maintaining  its 
strength  through  1700  GMT.  The  maximum  over  Mississippi 
loses  intensity  and  becomes  less  well  defined,  while  the 
maximum  originally  over  Georgia  moves  toward  the  southwest 
with  little  change  in  intensity.  An  additional  maximum 
develops  over  southern  Wisconsin  at  1600  GMT  (near  45)  and 
increases  in  strength  by  1700  GMT  (51). 

During  the  afternoon,  there  is  a  north-south  axis  of 
instability  along  the  Mississippi  River  Valley  at  2000 
GMT;  however,  this  changes  to  an  east-west  alignment  over 
Arkansas  and  Tennessee  at  2300  GMT  when  values  are 
greatest  (55).  Overall,  the  axis  of  instability  moves 
little  during  the  day  although  the  individual  maxima  show 
some  translation.  The  2000  GMT  analysis  appears  to  have 
the  best  agreement  with  convection  indicated  by  the 
satellite  imagery  (Fig.  6).  On  the  other  hand,  this  index 
also  fails  in  predicting  the  convection  over  eastern  Texas 
that  was  mentioned  earlier. 

K I  patterns  at  1300  GMT  (Fig.  24)  closely  resemble  the 
corresponding  LI  analysis  (Fig.  17).  Although  movements  of 
the  individual  KI  maxima  through  the  day  are  similar  to 
those  of  the  LI  representat i on ,  the  KI  fields  appear  to 
lose  definition  in  the  afternoon.  Furthermore,  they  do  not 
develop  the  maximum  over  Nebraska.  Thus,  KI  patterns  at 
2000  GMT  do  not  agree  as  well  with  VAS-derived  LI  (Fig. 
17),  and,  in  general,  they  do  not  appear  to  give  as  good  a 


representation  of  stability  as  the  LI  analyses.  This 
conclusion  for  21  July  agrees  with  that  of  Kitzmiller 
(1936)  who  found  the  KI  providing  less  information  than 
the  LI  during  the  afternoon. 

c.  Total  Totals  Index 

The  Total  Totals  Index  (TT)  is  commonly  used  by 
military  meteorologists.  It  is  quick  and  easy  to 
calculate,  requiring  information  from  only  two  levels.  The 
TT  is  given  by 

TT  =  (Tbho-  2T«„o)  +  DPooo •  (12) 

Radiosonde-derived  values  of  TT  normally  range  from  25  to 
50,  where  a  value  greater  than  40  indicates  the 
possibility  of  a  severe  weather  event.  Again,  VAS-derived 
versions  usually  are  larger,  ranging  from  30  to  65.  This 
index  also  can  be  modified  by  including  surface  data 
(Charba,  1979).  Since  the  modified  TT  gave  better  results 
than  the  original  version  (as  observed  with  the  KI),  only 
the  former  is  displayed  here. 

At  1300  GMT  (Fig.  25),  there  is  a  broad  north-south 
area  of  instability  stretching  from  Minnesota  to  the  Gulf 
of  Mexico.  By  1600  GMT,  an  axis  of  greater  instability 


devel ops 


over  central  South  Dakota  and  Nebraska  in 


response  to  the  VAS-anal yzed  500  mb  cold  trough  (Fig.  10). 
At  2000  GMT,  a  north-south  axis  of  instability  stretches 
from  northern  Illinois  to  western  Tennessee,  with  a 
maximum  o-f  68  over  central  Illinois.  Another  area  of  major 
instability  (68)  covers  eastern  Nebraska,  while  two  others 
(64)  are  located  over  eastern  Oklahoma  and  northwestern 
Florida.  At  2300  GMT,  the  instability  maximum  over 
Nebraska  becomes  the  dominant  feature. 

The  TT  analyses  (Fig.  25)  are  similar  to  those  of  KI 
(Fig.  24)  and  LI  (Fig.  17)  over  the  Midwest  and  Southeast. 
Also,  like  the  LI  analyses,  but  unlike  the  KI,  there  is  a 
strong  instability  region  near  Nebraska.  Kitzmiller  (1986) 
noted  that  the  TT  index  did  not  perform  as  well  as  the  LI; 
however,  there  is  no  clear  distinction  on  21  July. 


8.  Integrated  parameters 


The  previous  chapters  have  demonstrated  that  VAS  can 
reasonably  describe  horizontal  -fields  o-f  temperature  and 
dewpoint.  However,  since  weighting  functions  of  the 
sensor's  spectral  channels  are  so  broad  (Fig.  1),  vertical 
resolution  is  poor.  Several  authors,  e.g.,  Ilostek  al . 

(1986),  Zehr  (1985),  have  suggested  using  parameters  based 
on  data  from  deep  atmospheric  layers  rather  than  from 
discrete  levels.  Therefore,  this  chapter  evaluates  a 
variety  of  parameters  and  indexes  which  use  the  concept  of 
layer  integration.  Personal  communication  with  Dr.  James 
Purdom  of  the  National  Environmental  Satellite  Data 
Information  Service  (NESDIS)  was'  the  impetus  for  this 
aspect  of  the  research.  Dr.  Purdom  kindly  provided  a 
program  to  calculate  various  parameters,  and  it  was 
modified  by  the  author  to  derive  several  additional 
variables.  The  investigated  parameters  include: 

1.  Mean  mixing  ratio  in  lowest  kilometer, 

2.  Negative  energy  from  surface  to  LCL , 

3.  Pressure  of  LCL, 

A.  Height  of  wet-bulb  zero. 


.  Negative  energy  from  LCL  to  Level  of  Free  Convection 


(LFC )  , 

6.  Heating  energy  (HE)  required  for  surface  to  reach 
Convective  Temperature  (T<=), 

7.  Total  Negative  Energy  (NE)  , 

8.  Total  Positive  Buoyant  Energy  (PBE) , 

9.  Total  Positive  Energy  if  surface  temperature  equals 

-T- 

Tc  , 

10.  Maximum  Vertical  Velocity,  and 

11.  Precipitable  Water  (PW) . 

Zehr  (1985)  described  positive  buoyant  energy 
Convective  Available  Potential  Energy  (CAPE-  ,  .s 

vertically  integrated  temperature  difference  etwee 
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Fig.  26.  Skew  T-log  P  diagram  illustrating  bulk 
thermodynamic  parameters  computed  from  v'AS 
retrievals  (after  Zehr ,  1?S5). 

Area  A  represents  negative  area  from  sur¬ 
face  to  LCL. 

Area  8  represents  negative  area  from  LCL 
LFC. 

Area  C  represents  heating  energy  required 
for  surface  temperature  to  reach  Tc  . 

Area  D  represents  positive  buoyant  energy. 
Area  E  represents  available  positive  buoy¬ 
ant  energy  1 f  surface  reaches  Tc . 


wet-bulb  zero  is  self  explanatory.  The  heating  energy 
required  to  raise  the  surface  temperature  to  its 
convective  value  is  represented  by  area  C,  while  the  total 
negative  energy  is  merely  the  sum  of  parameters  #2  and  #5. 
Parameters  #8  (area  D)  and  #9  (area  E)  are  energies  that 
would  be  available  to  a  parcel  after  having  reached  its 
LFC  from  either  the  original  surface  temperature  (due  to 
forced  lifting),  or  from  its  convective  temperature, 
respectively.  The  maximum  vertical  velocity  is  obtained 
from  simple  parcel  considerations.  Finally,  precipitable 
water  is  integrated  water  vapor  content  throughout  the 
depth  of  the  sounding. 

The  eleven  parameters  were  calculated  at  1300,  1600, 
1700,  2000,  and  2300  GMT.  Although  each  provided 
reasonable  fields,  several  stood  out  from  the  others. 
Specifically,  positive  buoyant  energy,  total  negative 
energy,  heating  energy  required  to  reach  surface 
convective  temperature,  and  precipitable  water  analyses 
appeared  to  provide  the  most  useful  information.  As  a 
result,  they  will  be  described  in  the  following  pages. 

a.  Positive  buoyant  energy 

Fig.  27  displays  PBE  for  the  five  retrieval  times. 
Patterns  and  evolutions  of  these  analyses  are  consistent 
with  those  of  the  Lis  (Fig.  17).  For  example,  at  1300  GMT, 


major  energy  centers  are  located  over  eastern  Missouri, 


northern  Mississippi,  and  northern  Georgia.  During  the 
day,  PBE  increases  greatly  ( ^  3500  J/kg)  over  southern 
Alabama,  southern  Georgia,  and  northern  Florida,  while 
another  maximum  03000  J/kg)  develops  over  west  central 
Tennessee.  Satellite  imagery  (Fig.  6)  shows  the 
Southeastern  maxima  agreeing  well  with  ongoing  convection; 
yet  the  Tennessee  area  still  is  relatively  convection 
free.  Convection  does  take  place  shortly  thereafter. 

An  axis  of  maximum  PBE  extends  from  western  South 
Dakota  to  central  Nebraska  at  1700  GMT,  and  by  20CC  GMT 
it  becomes  a  strong,  well  defined  center  (3000  J/kg).  This 
is  similar  to  that  of  the  LI  analyses  (Fig.  17);  however, 
as  discussed  earlier,  no  convection  takes  place  until  much 
later.  Elsewhere,  a  minimum  (<  1500  J/kg)  develops  over 
eastern  Missouri,  persisting  through  the  late  afternoon. 
Finally,  relatively  large  values  over  eastern  Oklahoma  at 
1600  and  1700  GMT  and  the  minimum  axis  covering  portions 
of  central  Louisiana,  Mississippi,  and  Alabama  at  several 
times  both  agree  with  LI  features. 

Zehr  (1985)  noted  that  values  of  PBE  exceeded  750  J/kg 
in  areas  of  ongoing  convection  and  that  PBE  maxima 
outlined  preferred  areas  of  convection.  In  the  current 
case  however,  areas  >  750  J/kg  did  not  necessarily  have 
concurrent  convection,  e.g.,  most  of  the  Plains  States  had 
afternoon  values  greater  than  1500  J/Kg,  yet  skies 
remained  clear.  Obviously,  other  factors,  e.g.,  inadequate 


moisture  or  a  lack  of  -forcing  are  preventing  that  energy 
-from  being  realized.  Concerning  similarities  between  PBE 
and  LI,  this  is  expected  since  procedures  for  obtaining 
their  values  are  closely  related.  Specifically,  the  LI  is 
the  difference  between  the  lifted  parcel  temperature  and 
the  environmental  value  at  500  mb,  while  PBE  is  the 
integrated  temperature  difference  between  the  raised 
parcel  and  the  envi ronmental  sounding  at  all  levels  above 
the  LFC  (where  the  parcel  is  warmer).  Furthermore,  since 
VAS  data  only  are  available  at  the  mandatory  levels,  the 
500  mb  value  has  a  major  impact  on  the  overall  profile. 

b •  Total  negative  energy 

Total  negative  energy  describes  inhibition  to 
convective  activity.  The  negative  region  acts  as  a  hurdle 
which  the  sub-cloud  parcel  must  overcome  by  heating  or 
lifting  before  even  large  PBE  can  be  utilized.  Fields  of 
NE  are  presented  in  Fig.  28.  Whereas  PBE  (Fig.  27)  and  LI 
(Fig.  17)  analyses  compared  closely,  the  NE  patterns  are 
much  different.  In  this  case,  maximum  NE  represents  a 
stronger  "cap",  i.e.,  a  greater  amount  of  lifting  or 
heating  is  required  for  a  parcel  to  reach  its  LFC. 

At  1300  GMT  (Fig.  28)  ,  an  axis  of  minimum  NE  (  :  200 
J/kg)  extends  from  southern  Minnesota  to  eastern  Missouri, 


then  eastward  to  central  Kentucky,  and  finally  southward 


rvi 


0600  GMT. 


Overall,  there  is  a  consistent  relation  between  NE 
(Fig.  2B>  and  thunderstorm  activity  (Fig.  6).  That  is,  NE 
maxima  correspond  to  little  convection,  regardless  of  PBE 
strength,  and  rarely  does  deep  convection  occur  with  NE 
greater  than  250  J/kg.  This  alone  appears  to  be  a  valuable 
forecasting  tool. 

The  use  of  energy  parameters  for  predicting  convective 
activity  is  not  new.  Normand  (193S)  categorized 
atmospheric  soundings  into  three  classes:  no  latent 
instability,  psuedo-1 atent  instability,  and  real  latent 
instability.  When  PBE  was  zero,  the  sounding  was  said  to 
have  no  latent  instability.  Furthermore,  if  PBE  was 
greater  than  zero,  but  less  than  NE,  the  sounding  had 
"psuedo-1 atent "  instability.  On  the  other  hand,  if  the  PBE 
exceeded  NE,  then  the  sounding  had  real  latent 
instability.  Thus,  Normand  (193B)  attempted  to  combine  NE 
and  PBE  into  one  index  which  would  predict  convection; 
however,  his  theories  were  not  widely  accepted.  Others 
felt  that  combinations  were  not  relevant,  since  even  small 
negative  energy  can  prevent  a  large  positive  value  from 
producing  storms.  On  the  other  hand,  the  ability  to 
describe  these  parameters  at  the  smal 1  scale  (spatially 
and  temporally),  and  perhaps  more  importantly,  their 
changes,  demands  renewed  investigation  of  this  concept. 


Heating  to  reach  convective  temperature 


i 

I 

c  • 

There  are  several  mechanisms  which  allow  parcels  to 
overcome  their  negative  energy  barriers,  e.g.,  frontal 
lifting,  orographic  lifting,  and  heating.  Fig.  29  depicts 
the  heating  required  to  raise  the  surface  temperature  to 
the  convective  temperature,  i.e.,  the  amount  of  heat 
energy  which  would  raise  the  parcel's  LCL  to  its  LFC.  Like 
the  NE  analyses  (Fig.  23),  large  HE  indicates  that  a 

greater  amount  of  heating  is  nessesary  for  convection, 

therefore,  convection  based  on  heating  alone  would  not  be 
as  likely.  Conversely,  minimum  HE  indicates  a  preferred 
area  for  heating  induced  convection. 

HE  patterns  at  1300  GMT  (Fig.  29)  contain  a  minimum 
axis  (200  J/kg )  stretching  from  Indiana  southward  to 

northern  Florida,  while  a  maximum  covers  northeastern 

Texas  (700  J/kg).  The  Northeast  is  a  region  of  large  HE. 
HE  decreases  over  the  entire  area  by  1600  GMT  (Fig.  29) , 
and  this  is  expected  due  to  diurnal  heating.  Although  the 
minimum  axis  does  not  appear  to  change  position,  a  maximum 
axis  moves  out  of  the  Dakotas  into  Nebraska  and  Iowa.  By 
1700  GMT,  the  north-south  minimum  remains 

quasi -stati onary;  however,  a  tongue  of  minimum  values 
develops  across  northern  Arkansas  and  southern  Missouri 
(150  J/kg).  This  particular  axis  seems  to  match  the 
outflow  boundary  from  the  morning  convection  over  central 


Missouri  that  was  suggested  by  satellite  imagery  (Fig.  6). 


The  maximum  axis  appears  to  strengthen,  wi th  a  center 
developing  over  Oklahoma  (650  J/kg). 

During  the  late  afternoon  hours  (2000  and  2300  GMT, 
Fig.  29),  the  major  minimum  axis  continues  to  cover  the 
Southeast;  however,  the  northern  portion  appears  to  move 
slightly  westward  into  northern  Illinois.  Most  of  the 
maxima  show  only  minor  fluctuations;  however,  values  over 
Nebraska  fall  507.,  comparable  to  results  of  NE  (Fig.  23). 

There  does  not  appear  to  be  a  well  defined  "critical'1 
HE  value  which  delineates  the  presence  or  absence  of 
convection.  For  example,  thunder storms  form  over  Arkansas 
where  HE  reaches  375  J/kg  (recall  the  outflow  boundaries 
in  Fig.  6),  and  they  develop  over  eastern  Texas  where  HE 
exceeds  600  J/kg.  This  may  suggest  that  HE  predicts  areas 
of  preferred  convective  development,  but  will  not  indicate 
areas  where  convective  growth  is  inhibited. 

d .  Precipitable  water 

Water  vapor  content  is  an  ingredient  to  many  of  the 
stability  parameters  described  here  and  in  Chapter  7,  and 
it  is  valuable  to  consider  it  separately.  Precipitable 
water  represents  the  integrated  vapor  of  an  atmospheric 
column.  Fig.  30  displays  values  for  the  total  surface  to 
300  mb  1  aver .  Values  range  from  2.5  to  6.5  cm,  which  is 
typical  of  RfiOB - der i ved  values  during  summertime.  Patterns 


and  evolutions  are  similar  to  those  of  S50  mb  dewpoint 
(Fig.  9).  An  axis  o-f  maximum  PW  (>  4.0  cm)  extends  from 
southern  Wisconsin  to  the  Gulf  Coast  at  1300  GMT  (Fig. 
30),  shifting  slowly  to  the  east  with  time.  Values 
reaching  6.0  cm  develop  over  Mississippi  at  1600  GMT  and 
move  north  during  the  afternoon,  eventually  covering 
southern  Illinois  where  the  intense  convection  develops. 
Values  are  less  than  2.0  cm  over  the  Northeast  and  the 
eastern  Rockies,  the  latter  mostly  due  to  the  reduced 
thickness  of  the  column. 

To  explore  vertical  characteri sties  of  the  moisture, 
PW  was  calculated  over  three  component  columns:  Surface  - 
700  mb,  700  -  500  mb,  and  500  -  300  mb.  Results  for  2000 
GMT  are  shown  in  Fig.  31.  As  one  might  expect,  most  of  the 
moisture  is  found  in  the  lowest  layer.  The  middle  layer 
appears  to  contain  about  257.  of  the  total  amount  while  the 
upper  layer  contains  less  than  5  /..  There  appears  to  be  a 
slight  eastward  shift  in  the  PW  axis  with  altitude,  with 
the  upper  level  position  extending  approximately  100  km 
east  of  the  low  level  axis  over  Illinois.  It  is  possible 
that  cirrus  blowoff  from  the  afternoon  convection  over 
western  Illinois  already  had  advected  to  the  east, 
contributing  to  the  vertical  slant  of  the  moisture  axis. 
The  low  levels  over  eastern  Missouri  contain  a  minimum  of 
vapor  due  to  the  outflow  boundary;  however,  no  major  areas 
of  mid  level  dryness  are  evident  in  PW  or  in  the 


individual  soundings  (Figs.  18-22) .  Consequently, 
convective  instability  is  not  a  -factor  in  producing  the 
storm  outbreak.  Overall,  PW  fields  differ  only  slightly 
from  those  of  850  mb  dewpoint,  and  they  provide  little 
additional  information. 


9.  Summary  and  conclusions 


A  "new"  physical  retrieval  algorithm  which  converts 
atmospheric  radiances  measured  -from  a  space  platform  into 
vertical  profiles  of  temperature  and  moisture  has  been 
evaluated  for  a  "typical"  summertime  day  in  July  1982.  The 
area  of  interest  was  the  midwestern  United  States  where 
data  coverage  was  excellent.  Sounding  products  were 
available  at  intervals  of  2-3  h  for  a  total  of  5  times. 
After  a  brief  review  of  the  two  different  retrieval 
algorithms,  RMS  errors  for  the  new  simultaneous  scheme 
were  compared  to  those  from  an  earlier  version,  using 
radiosonde  reports  as  a  standard.  VAS-deri ved  fields  of 
temperature  and  dewpoint  were  evaluated  for  both  spatial 
and  temporal  continuity.  Next,  numerous  stability  indices 
calculated  from  the  satellite-derived  profiles  were 
compared  to  visible  satellite  imagery.  Their  ability  to 
define  small  scale  stability  features  and  information  not 
obtainable  from  conventional  sources  was  examined 
ex  tensi vel y . 

Overall,  soundings  from  the  new  retrieval  algorithm 
were  favorable;  however,  there  also  were  some  limitations. 
Compared  to  results  for  the  older  version,  there  was 
marked  improvement  in  moisture  analyses  which  exhibited 


good  continuity  -from  one  time  to  another.  On  the  other 


hand,  while  temperature  patterns  corresponded  to  observed 
RAOB  analyses,  mean  errors  were  somewhat  greater  although 
RMS  errors  changed  little.  The  newer  VAS-derived  soundings 
still  exhibited  errors  considerably  greater  than  their 
RAOB  counterparts ,  and  unrealistic  diurnal  trends 
continued  to  occur.  The  theoretical  aspects  of  how  changes 
in  the  algorithm  have  altered  its  overall  performance  were 
not  the  subject  of  this  research.  However,  it  appears  that 
the  mare  realistic  "first  guess"  water  vapor  profile  is  a 
significant  contributor  to  the  noted  improvement  in 
moisture  analysis.  The  problem  of  biased  radiances  from 
the  radiometer  underscores  the  need  to  consider  retrieval 
variations,  and  not  actual  magnitudes. 

Nan-i ntegrated  and  integrated  stability  indices  were 
calculated  because  of  the  improvement  in  moisture 
definition.  The  Lifted  Index  appears  to  be  a  promising 
VAS-derived  convection  predictor.  LI  anal yses  maintained 
reasonable  continuity,  and  provided  useful  information 
about  several  locations  which  was  not  available  from 
either  RAOB  or  surface  observati ons.  In  many  cases,  the 
VAS-derived  LI  patterns  were  superior  to  those  from  RAOBs; 
however,  the  satellite's  values  were  considerably  more 
negative.  Apparently,  better  low-level  moisture  analyses 
were  a  major  contributor  to  this  success. 

Total  Totals  and  K  Indexes  appeared  to  give  a 


convection  forecaster  reasonable  information,  but  patterns 
were  not  as  well  defined  as  those  of  LI.  Also,  most 
features  in  the  KI  and  TT  analyses  were  present  in  the  LI 
fields.  In  other  words,  KI  and  TT  often  detected  the  same 
major  stability  centers  as  did  the  LI;  however,  the  LI 
analyses  better  defined  smaller  scale  features,  e.g., 
minor  axes  of  stablility. 

Eleven  integrated  parameters  were  calculated;  several 
provided  valuable  information.  Speci f i cal  1 y ,  fields  of 
positive  and  negative  buoyant  energy  appeared  to  isolate 
areas  which  later  contained  strong  convection,  or  on  the 
other  hand,  remained  convection  free.  As  expected,  there 
generally  was  good  agreement  with  the  non-i ntegrated  LIs. 
Heating  energy  required  to  raise  the  surface  to  its 
convective  temperature  provided  some  insight  into  the 
"timing"  (onset)  of  convection;  however,  these  fields 
generally  were  si  mi  liar  to  those  of  NE.  Finally,  VAS 
precipitable  water  analyses  closely  matched  low-level 
moisture  fields  depicted  by  VAS,  and  values  were 
consistent  with  RAOB  analyses. 

It  may  seem  that  the  author  is  overly  optimistic  about 
the  performance  of  VAS-derived  parameters.  However, 
having  spent  several  summers  forecasting  convection  over 
the  United  States  at  AFGWC  in  conjunction  with  the  NSSFC, 
it  quickly  became  clear  that  in  the  absence  of  stronger 
springtime  dynamics,  a  good  summertime  convective  forecast 
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